Holographic Super fluids and Superconductors 

in Dilaton-Gravity 



Alberto Salvio 

Scuola Normale Superiore and INFN, Piazza dei Cavalieri 7, 56126 Pisa, Italy 

(N ■ 
t— i ' 

O ' Abstract 

(N ■ 

We investigate holographic models of superfluids and superconductors in which the gravita- 
. tional theory includes a dilatonic field. Dilaton extensions are interesting as they allow us 

C/3 " to obtain a better description of low temperature condensed matter systems. We focus on 
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,— 1| infrared and is therefore compatible with the normal phase being insulating (or a solid in 

■ the fluid mechanical interpretation); indeed we find that this is the case at low temperatures 

JL . and if one appropriately chooses the parameters of the model. Not only the superfluid phase 
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1 Introduction and results 



The anti de Sitter/conformal field theory (AdS/CFT) correspondence is a powerful method to 
investigate strongly coupled theories. An interesting application of such duality is the study of 
systems at finite temperature and density of some conserved quantum number. This is the typical 
situation in condensed matter physics and in the last few years a considerable amount of work has 
been dedicated to the development of holographic methods for condensed matter systems (reviews 
of which can be found in Refs. [1, 2, 3, 4, 5, 6, 7]). 

In the simplest realizations one deals with Einstein's gravitational theory with a negative 
cosmological constant coupled to a gauge field (dual to the conserved current). While this setup 
captures some interesting properties of realistic systems, the low temperature behavior has some 
disturbing features. A classic example, which is perhaps unphysical rather than disturbing, is 
the fact that the simplest charged black hole solution predicts a nonvanishing entropy in the zero 
temperature limit. This situation can be improved [8, 9, 10] by adding a real scalar field (a 
dilaton), which interacts with the gauge field by modifying its gauge coupling g, 

i _L %M m 

9 2 9 2 (<P) 9 2 ' 1 ] 

where 7^ 0, and self-interacts through a potential. This allows for zero entropy at zero tem- 
perature. Moreover, real scalars of this type are typically present in the low energy limit of string 
theories and their inclusion may be good if eventually one has to find a stringy embedding. 

In this work we initiate the study of holographic superfluids and superconductors in dilaton 
extensions of Einstein's gravity. Like in the absence of the dilaton, the minimal description of 
these systems requires a conserved U(l) current and a charged condensing operator, which are 
respectively dual to an Abelian gauge field and a charged scalar living on an asymptotically AdS 
space [11, 12, 13, 14] (for a review see [1, 15, 16, 6]). Then the dilaton can couple directly not only 
to the gauge field but also to the charged scalar 1 . In the case of superconducting materials the 
U(l) symmetry is gauged while for superfluids is global. In holography a superfluid corresponds 
to a Dirichlet condition for the gauge field at the AdS boundary, while a true superconductor 
can be obtained by imposing a boundary condition of the Neumann type [18]. Superfluids can be 
considered as superconductors in the limit of nondynamical (electromagnetic) gauge fields. For 
the sake of definiteness, we adopt the terminology used in the literature of superconductivity, 
rather than that of superfluidity. 

One of the main motivations for holographic superconductivity is the hope to shed light on 
high-T superconductors. Experiments indeed imply that these materials cannot be entirely de- 
scribed by the weakly coupled theory of Bardeen, Cooper, and Schrieffer (BCS), but they could 
involve strong coupling. The cuprate high-T superconductors can be obtained by doping a Mott 
insulator and in the (T versus "dopant concentration" ) phase diagram one can observe an insula- 
tor/superconductor transition in the low temperature region [19]; therefore the normal (nonsuper- 
fluid) phase is expected to be gapped. A transition of this sort has been realized in holography by 
considering a theory with a compactified extra dimension [20, 21] (see also [22] for the realization 
of a true superconductor with a dynamical gauge field). The size of this extra dimension breaks 



1 Our setup is different from that in Rcf. [17] because the dilaton there is defined as the modulus of the charged 
scalar rather than a real scalar, as it is usually done. 
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scale invariance and at low temperatures leads to a gap of charged excitations and thus to an 
insulator normal phase. One expects to be able to obtain an insulator/sperconductor transition 
by introducing a dilaton, even in the absence of additional dimensions, because a running dilaton 
also introduces a scale. We will find that this is indeed the case at low enough temperatures and 
for some values of the parameters. 

Let us describe in more detail our main results and the organization of the paper. In section 2 
we will define our Einstein-dilaton model coupled to the U(l) gauge field and the charged scalar 
and recall the AdS/CFT dictionary. Here and in the rest of the paper we focus on 2+1 dimensional 
CFTs. In the same section we will also review the (recently found [23]) most general static and 
asymptotically AdS planar black hole of the dilaton-gravity system. 

Once we fix the parameters of the model, there is only one solution with active dilaton (dilaton- 
BH), but the fact that the space is asymptotically AdS implies that the Schwarzschild black hole 
(S-BH), which has <fi = 0, is also a solution. Therefore, in section 3, we compare their free 
energies as a function of T to determine which one is energetically favorable. We find that the 
S-BH is favorable at high temperatures, while the dilaton-BH dominates in the low temperature 
region. In the latter phase the dilaton runs logarithmically with the energy and grows from the 
ultraviolet (UV) to the infrared (IR). We then study the AC and DC conductivity by considering 
an electromagnetic (EM) perturbation on top of these geometries. The presence of a black hole 
horizon implies that the system is never perfectly insulating, but for some choices of Za{4>) the DC 
conductivity is suppressed especially at low temperatures, resembling an insulator. This occurs 
when the running coupling constant g 2 (4>) becomes large in the IR and, not surprisingly, links the 
insulating phase with IR strong coupling. 

In section 4 we find superfluid phase transitions by solving the gauge field and charged scalar 
equations on the fixed dilaton-BH geometry. Like for the S-BH [12], for each value of the tem- 
perature and the other parameters there is a critical temperature T c below which the charged 
operator condenses breaking the U(l) symmetry. We also study the conductivity in the broken 
phase and obtain that both the gap of charged excitations and the superfluid density are modified 
by changing the couplings of the dilaton, in particular the running constant, Eq. (1). 

Some of the most famous properties of superconductors occur in the presence of external 
magnetic fields; therefore we move on and study the magnetic response in the dilaton-BH phase. 
The same type of analysis in the model without the dilaton has been performed in [14, 24, 18, 22]. 
While for superfluids the external and total magnetic fields coincide, they differ substantially in 
true superconductors. In section 5 we study the superconductor case. We directly observe the 
Meissner effect and find vortex solutions, the penetration depth and the critical magnetic fields H c i 
and H C 2. Like in the case with no dilaton [18, 22], we always find H c \ < H C 2 (the superconductor is 
of Type II) and thus this seems to be a rather general prediction of holography. Quite interestingly, 
the known high-T superconductors are of this type. Crucial differences with respect to the S-BH 
phase [18] emerge in the low temperature region, for some choices of Za{4>) (which, interestingly, 
lead to an insulating normal phase). In these cases the ratio H c2 /H c i remains finite in the limit of 
small T/T c unlike the S-BH case [18] and also the penetration depth shows an analogue difference. 
Finally in section 6 we study superfluids in the presence of external rotations and find superfluid 
vortices in the dilaton-BH phase; here we find results that are qualitatively similar to those for 
the S-BH [25, 26, 18], although the dilaton generically modify the vortex shape. 

Section 7 contains some outlook of our work. 
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2 The model and the AdS/CFT correspondence 



Our model is defined by the usual Einstein-dilaton action coupled to a U(l) gauge field A a and a 
charged scalar 

s = j^{^[n-^?-vm -^. 8 -^W}, , 2) 

where TZ is the Ricci scalar, Gn the gravitational Newton constant, a real scalar (the dilaton) 
and we introduced T a p = d a Ap — dpA a and D a = d a — iA a . A negative cosmological constant A 
is included in the dilaton potential V(4>), and the AdS radius is defined by A = —6/L 2 . We define 
such that V(0) = A. Also the dilaton couples to A a and \1> through two generic functions Za(4>) 
and Z^((f>) and we choose the normalization of the matter fields in a way that Za{0) = Z^(0) = 1. 
There are no special requirements for the Zs at this level, besides the fact that they should be 
regular and nonvanishing for any <fi in order for the semiclassical approximation to be valid. For 
simplicity, we have not added any potential for the charged scalar. 

The most general static asymptotically AdS planar black hole with two-dimensional rotation 
and translation invariance has the form 

ds 2 = W{z) (-f(z)dt 2 + dy 2 + ^ , <P = cf>(z), (3) 

where t represents time, dy 2 = 5ijdy l dy : > and z is the holographic coordinate, which we define 
such that the AdS-boundary is located at z = 0: for z close to zero, W, f and <fi approximate 
to W{z) = L 2 /z 2 , f(z) = 1 and </>(z) = 0. This allows us to use the standard AdS/CFT 
dictionary, which relates the property of a gravitational theory with those of a CFT. In particular 
the boundary values of \& and A^ are the sources of a charged operator O and the current J M in 
the CFT side. Therefore, if one solves the classical field equations with boundary conditions 

#|, =0 = tf , \\z=o = a„, (4) 

the Green's function for O and are given by functional derivatives of the on-shell action with 
respect to \l/o and a M . For example the vacuum expectation values are given by 

(J M ) = 4^l*=o , (O) = -^D z V\ z=0 . (5) 
9 9 z 

We will work in the limit Gn — > such that the effect of the gauge field and the charged scalar 
on the metric and the dilaton can be neglected. In this limit the most general asymptotically AdS 
black hole of the form given in Eq. (3) is [23] 



(-V(tA + (i + ^ 



z, I I 4 — v 2 



{l + z/Lf (1 + z/L)-" 



V + l l — v 



(6) 
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The black hole horizon Zh, defined as usual by f(zh) = 0, is a function 2 of v and z s . Since 
we want to study the theory at finite temperature, we perform the Euclidean continuation to 
imaginary time it with period (3 = 1/T. The temperature can then be computed with the formula 
T = \ f f (z h )\/4:7i, where a prime denotes a derivative with respect to z. 

When v — 1 we have = and we recover the S-BH: W{z) = L 2 /z 2 ,f(z) = 1 — (z/z s ) 3 . 
In this case z s = z^ and we have Zh = 3/(AnT). If instead v ^ 1 the dilaton is active and its 
backreaction on the metric is nontrivial. This results in a displacement of the horizon and the 
temperature is no longer given by 3/(47rz s ), but by a more general (z/-dependent) expression. It 
is important to notice that the S-BH is a solution whenever V has a negative stationary point. In 
order to have the dilaton-BH, v ^ 1, one has to consider instead the following family of potentials 
[23] 



Vffl = £ v&r** (7) 



where Vi and Si are given by 



z l V . = o l - u i 2 }- u2 '1 ■ ,^ 2 -(^ 2 -4)zf/L a 



L 2 + i/' i/ 2 -4' v ' z/2(z/-2) 

zz+1 3^-(^-4)^/L3 n 3, 2 -(, 2 -4)^3 /L 3 
3 2~^' 4(Z/ " 1} z/2 (zy2 _ 4) . 7^7^ . (8) 



and 



z/ 2 



=Fy— J— *i = f + 1, i, i-v, -i, -v-i (9) 

respectively, where the + signs correspond to > and <ft < 0; the potential is therefore symmetric 
under <ft — > —<f) and from now on we restrict to > 0. Since the black hole is asymptotically AdS, 
this potential has a stationary point for = such that V(0) = A. An immediate but important 
consequence is that one has two black holes for any v ^ 1, the S-BH and the one with active 
dilaton. For small enough one obtains V((f>) ~ A — -^-0 2 . Notice that the squared mass of is 
always above the Breitenlohner-Freedman bound [27] and therefore = is never unstable, but 
either metastable or stable. 



3 The Einstein-dilaton theory 

In this section we focus on the normal phase, (O) = 0. This, among other things, will help us to 
understand what kind of superfluid transitions are those described in the rest of the paper. 

3.1 Phases 

Let us compare the free energies of the S-BH and the dilaton-BH to determine which one is 
energetically favorable, as a function of T. The free energy F can be computed by evaluating the 
euclidean action Se at the solution of interest, 

F = (3S E - (10) 

2 By dimensional analysis z/,/L is a function of v and z s /L only. 
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According to this definition it is clear that the free energy of the S-BH is divergent. One reason 
is that the solution does not depend on the planar coordinates y and therefore F goes to infinity 
as the volume of the plane V = J d 2 y — > oo. This is an infrared divergence and can be avoided 
by considering the free energy density F/V. There is, however, another source of singularity, due 
to the fact that the Lagrangian evaluated at the S-BH goes as 1/z 4 for z small enough and the 
integral includes the point z = 0; however this is an ultraviolet divergence and as such it should 
cancel in the difference Fj — F s , where F^ and F s are the free energy densities of the dilaton- 
and S-BH respectively. Indeed the dilaton-BH is asymptotically AdS and therefore its ultraviolet 
behavior should coincide with that of the S-BH. 

We plot the free energy difference in Fig. 1. We find that if we lower the temperature below a 
critical value T c there is a phase transition between the S-BH and the dilaton-BH. Since L is the 
only dimensionful parameter, % oc 1/L, but the actual value of % depends on the dimensionless 
parameter v as well. 
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Figure 1: The free energy density (up to the gravitational constant 1/16ttGn ) of the dilaton-BH minus 
that of the S-BH as a function of the temperature. 

We observe that the entropy, —8tF, approximates to zero at small temperatures. 
3.2 Conductivity 

Both the S-BH and dilaton-BH phases correspond to systems with nonvanishing conductivity in 
the CFT side. This property can be viewed essentially as a consequence of the black hole horizon 
(see below). However, as we will see, the DC conductivity can be suppressed for some values of 
the parameters, resembling an insulator. 

To compute the conductivity let us consider, on top of these geometries, a small time-dependent 
perturbation along a spatial coordinate x, 

A x (t,z)=A(z)e l ^- t \ (11) 

which is induced by a small electromagnetic field at the AdS-boundary, a x (t). Here A and p are 
real functions of z. The system responds creating a current which is linear in a x : (J x ) = o~E x , 
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where a is the complex conductivity and E x = —dta x the electric field. Using the AdS/CFT 
dictionary, the first equation in (5), we have 

(12) 

Since this is a linear response problem the conductivity can be computed by solving the linearized 
Maxwell equation 

d z (fZ A (cf>)d z A x ) + W 2 ^A, = (13) 

with appropriate boundary conditions. At the AdS-boundary the condition is, as we said, A^l^o = 
a x , while at the black hole horizon we impose that the EM wave is going towards the horizon 
(ingoing boundary condition) 

p'(z)f(z) 41 at z = z h , (14) 

as required by regularity [1]. 

Inserting (11) into (13) we obtain an equation whose imaginary part can be solved analytically, 
to obtain 

1/(0) (A(0)V 
P[Z) ~ Z A (4>(z))f(z){A(z)J ■ (15) 

A consequence is that p'(0) cannot be zero without violating Eq. (14). Therefore the amplitude 
of the EM wave, A{z)e luip ^ z \ acquires an imaginary part which generates a nonvanishing Re [a]. In 
the fluid mechanical interpretation Re[cr] ^ corresponds to a fluid, as opposed to a solid which 
has vanishing Refer], as discussed in [22]. 
The real part of Eq. (13) is instead 



d z (fz A d z s) + fz A {d z sf + ( l - %^ ) = o, (16) 



-4s 

/ Z A 

where 



to be solved with the regularity boundary conditions 

1 , 2u?Z' A jZ A 

s = --\nZ A , s =-___, a t z = z h . (17) 

Once a solution to this equation is found we can compute the real and imaginary part of the 
conductivity through 

e ~2s(o) g/(Q\ 
Re[er] = —J-, lm[a] = — p, (18) 

where we used Eqs. (12) and (15). 

A particularly simple case is when there is no direct coupling between the dilaton and the 
gauge field, Z A = 1, which corresponds to s = and thus to a frequency-independent Refer] and 
to a vanishing Im[<r], both for the S-BH and the dilaton-BH. In this case the dilaton does not have 
an impact on cr; it coincides with that found in [12]. 
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When Z4 depends on the situation is different. A simple way to see it is to take the zero 
frequency limit. Then Eqs. (16) and (17) tell us s(z) = —^hiZ^\ z=Zh for every z, which leads to 

lim Refer] = \z A \ z=Zh , limlrn[cr] = 0. (19) 

Here we see that the DC conductivity can be suppressed or enhanced depending on the function 
Za(<P). For the exponential form Za{4>) = e 7< ^, the more positive 7 is the bigger the DC conduc- 
tivity, while a large negative value of 7 corresponds to an approximate insulating behavior. This 
effect is even stronger at low temperatures where the dilaton is logarithmically large. We have 

lim Re[a] ~ T^v^ 2 " 1 )/ 2 at low temperatures, (20) 

u— >0 

where we have remembered that we restrict to the positive sign choice in (6), without loss of 
generality. The conductivity can depend on T because there is an additional scale in the dilaton 
phase. It is interesting to notice that for the functions Za{4>) such that the DC conductivity is small 
the phase transition described in Section 3.1 is between a conductor at high temperatures and a 
(non ideal) insulator at low temperatures, while in the other cases we observe a transition between 
two metallic phases; in particular if we set 7 = ^2/{y 2 — 1) in Eq. (20) we get a low temperature 
resistivity that goes linearly in T, like for many unconventional (not BCS) superconductors. 

In the high frequency limit, instead, Eqs. (16) and (17) tell us that s — > — | lnZ^ and so Eqs. 
(18) imply that a converges to the constant value we have in the S-BH phase: Re[cr] = 1/g 2 and 
Im[<r] = 0. At high enough energies these systems respond in a universal way. 

In Fig. 2 we give a as a function of 00. A direct coupling between the dilaton and the bulk 
gauge field gives us a nontrivial frequency dependence of both Re[cr] and Im[cr], unlike the S-BH 
phase (and the case Za{4>) = 1)- As suggested by the first equation in (19), Re[cr(u;)] is large 
(small) when Z A {4>) is large (small) at the horizon. We observe that the same is true for Im[er]. 
However, while for lm[a] the effect of the dilaton coupling function is strongest at intermediate 
values of u (as implied by the analytic discussion above), for Refer] we observe the maximum effect 
at low frequencies. 
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Figure 3: The condensate (O) as a function of fiL. We set v = 1.2, Z A {<t>) = e 7 ^ and = e act> . We 
choose T ~ 0.15/L while T ~ 0.010/L in the inset. 

4 Superfluid phase transitions 

In this section we obtain superfluid phase transitions by introducing a chemical potential. At low 
temperatures, for the functions Za(4>) such that the normal phase resembles an insulator, these 
are new insulator/superconductor transitions. 

4.1 Homogeneous superconducting phase 

In order to obtain superconducting solutions the simplest ansatz is 

* = A = A (z), (21) 

where is a real function, and the other components of A a are set equal to zero. A nonvanishing 
value of Ao is needed in order to have a regular solution for the condensate profile ip(z) (see below). 
Therefore we consider the following AdS-boundary conditions (see Eq. (4)) 

* = 0, a = fi. (22) 

The first condition ensures that the U(l) breaking is spontaneous, while the chemical potential \i 
is needed to have A ^ 0. The regularity at the black hole horizon indeed requires A (zh) = (in 
order for Aq to have finite norm) and therefore ao = would imply Aq = everywhere. 
The field equations for the ansatz in (21) are 

^d z (WfZ4<p)d z ^) + ^A& = 0, d z {Z A ^)d z A ) - 2 ^ (0 V A) = 0. (23) 

we can now see that A is needed: setting A = would imply either ip = or ip(z) oc 
J Q Z dz' I (W f Z^p) , which is singular at the black hole horizon. The requirement of regularity gives 
the following boundary conditions at Zh 

d z ^ = 0, A = 0. (24) 
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The condition Ao = ensures that A a has finite norm [12]. The second equation in (23) and the 
regularity of the solution implies that Aq should go to zero at least as fast as Zh — z and then 
the second term in the first equation of (23) should vanish at z = z^. The remainder of that first 
equation is a term proportional to /, which goes to zero for regular solutions, and one proportional 
to /' that can be easily simplified to obtain the first condition in Eq. (24). 

We have solved numerically the field equations in (23) with the UV and IR boundary conditions, 
in (22) and (24) respectively. In Fig. 3 we give the condensate as a function of the chemical 
potential. 

Like in [12] there is a critical chemical potential /i c above which O condenses, but here fi c /T 
depends on TL, v and the functional form of Z A and Z^. Inverting this relation one has the 
critical temperature of the superfluid phase transition, T c . We have that T c /li is a function of liL, 
v, Za and Z$. By computing explicitly the free energy, Eq. (10), one can see that, when it exists, 
the superconducting solution is energetically favorable with respect to the normal one. 

We observe that the condensate in Fig. 3 seems to diverge as iiL — > oo. Since increasing jxL 
for TL fixed corresponds to exploring the region of small T/T c , this suggests that L 3 (0) becomes 
bigger and bigger by decreasing T/T c . We notice that neglecting the backreaction of the charged 
scalar and the gauge field on the metric is still consistent in the strict Gn — > limit, which we 
are using throughout this paper 3 . 



4.2 Conductivity in the broken phase 

We now compute the conductivity in the spontaneously broken phase, (O) ^ 0. Most of the 
equations given in Section 3.2 are still valid here, but some of them have to be modified because 
the order parameter is nonvanishing: the linearized Maxwell equation, Eq. (13), now becomes 

d z (fZ A (<f>)d z A x ) + w 2 ^4 - ^Z^ 2 A X = 0. (25) 

In particular the equation for the modulus of the EM wave, parametrized by s, is modified to 

d z (fZ A d z s) + fZ A (d z s) 2 + u 2 ^ 1 - — - - — Z^ 2 = (26) 



/ V Z\) L? 
and has to be solved with the regularity boundary conditions 

1, 7 / V'WZ^ 2 /L 2 - 2uj 2 Z' a 
s = --\nZ A , s= v , at z = z h . (27) 



2 " Z A (Acu 2 + / 

At very high frequencies, Eqs. (26) and (27) indicate that the conductivity should converge to 
the value in the normal S-BH phase, regardless of the presence of the condensate. This corresponds 
to the physical fact that at energies much bigger than the condensation scale the system behaves 
as there were no condensate at all. For lower values of u we provide a in Fig. 4. The imaginary 
part diverges like l/uasu-^0, which corresponds, through the Kramers-Kronig relation 

Im[aM] = -hp r du' Re[ ^' )] , (28) 



3 We plan to go beyond the Gn — > limit in a future publication. 
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Figure 4: The conductivity as a function of the frequency in the broken phase, [i > [i c , (solid blue lines) 
and in the unbroken phase (dashed purple line). Re[o~] has also a term oc 5(u). We set T ~ 0.010/L, 
v = 1.2, Za{4>) = c 1 ^ and Z^ = 1. The upper plots have 7 = 0.4 (for which fi c ~ 0.37/L) while the lower 
ones correspond to 7 = —1 (fi c ~ 0.26/ L). 

to a delta function in the real part, Re[cr(u;)] ~ irn s 5(uj). This means that the system is supercon- 
ducting, or a superfluid in the fluid mechanical case. The quantity n s can be identified with the 
superfluid density. We also see that the presence of the condensate leads to a superconducting gap 
in the plot of Re[<r] versus u, like in the S-BH phase [12]. In the dilaton phase, however, the gap 
depends on Za(4>)- From the Im[er(u;)] plot we notice that the superfluid density is also changed 
by such function of the dilaton. 

5 Holographic superconductors in dilaton-gravity 

As we mentioned in the introduction, the distinctive property of superconductors is that the U(l) 
symmetry is gauged rather than global, as opposed to the superfluid case; all we have discussed 
so far can be applied equally well to superconductors and superfluids because the dynamics of 
the gauge field was not important 4 . The purpose of this section is to consider field configurations 

4 This statement refers to physical phenomena described by a„ and the condensate (a U(l) breaking field) only: 
when we say that the gauge field is (is not) dynamical we mean that the backreaction of the condensate on the 
gauge field has (has not) been taken into account, but we do not state anything about the interaction between a M 
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(mostly vortices) where the dynamics of the gauge field is instead crucial. When the gravitational 
theory is described by Einstein's gravity, this has been done in [18]; here we will consider the case 
of dilaton-gravity. 



5.1 Dynamical gauge fields in AdS/CFT 

In the standard AdS/CFT correspondence the boundary gauge field is considered as an external 
source for the U(l) current and as such is a nondynamical field, introduced through a Dirichlet 
boundary condition (the second equation in (4)). However, it is possible to promote to a 
dynamical field by substituting this Dirichlet boundary condition with one of the Neumann type, 
which for a 2+1 dimensional CFT is 



1 

? 



;.//( 



°ext 



0. 



(29) 



where J^ xt is an external nondynamical current and we have added a bare boundary kinetic term 



for generality [18]. This condition requires to add the following boundary terms to 



the action in (2) 



d 4 x 



(30) 



In the 2+1 dimensional case, which we are considering, the bare kinetic term in (29) is not 
necessary [28] and one can preserve conformal symmetry in the UV. In the absence of such term 
the dynamical a M can be considered as an emergent gauge field [18]. 

This method to introduce a dynamical gauge field has been originally applied to the holographic 
superconductor model [12] (including diverse dimensions [29]) in [18] and we refer to this paper for 
further details. It has been subsequently applied to holographic p-wave superconductors [30]- [33] 
in [34, 35]. It can be applied without any modification to the dilaton-gravity model of holographic 
superconductors presented here because the BH solutions we consider are all asymptotically AdS. 



5.2 Vortex phase 

Let us now look for genuine superconductor vortices with this recipe in mind. The simplest ansatz 
for the bulk fields is [36, 25, 37] 

= ^{z,r)e m \ A = A ( Zj r), A ip = A^(z 1 r) (31) 

and the other components of A a set equal to zero. In Eq. (31) r and ip are the polar coordinates, 
dy 2 = dr 2 +r 2 dtp 2 , restricted to < r < oo and < tp < 2tt, and n is an integer. The corresponding 

and the condensate constituents (whose operators arc not explicitly introduced here). 
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field equations are 



^d z (wfz 4 ,(m^ + ^W^) + (f - {A \ 2 71)2 



^(/Z40)d^) + Z A (<^ = 0, (32) 

d z (Z A (mA ) + ^dArd r A )- 2 -^j^^A = 0. 

Requiring regularity we obtain the following boundary conditions at z = Zh (with a reasoning 
similar to that below Eq. (24)) 

f{z h )d^ + -d r {rdri>) - ^ ~ = 0, 

f'(z h )d z A v + rd r (^<W) - ™Z A (t) ^ {A * ~ n) = °' (33) 



A = 0, 



and at r = 



<9 r A = A^ = 0, 
<9 r V> = for n = 0, V 7 = for n ^ . (34) 

For r — > oo we impose instead the physical condition that the fields should approach the homo- 
geneous superconducting configuration: 

d r ip = , <9 r A = , A v = n. (35) 

Finally let us consider the boundary conditions at z — 0. Requiring the U(l) to be spontaneously 
broken we impose (22), while demanding the magnetic field B = d r a^jr to be dynamical we have, 
from Eq. (29), 

1 „ fl 
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+ ~^rd r -d r A v 
2=0 el \r 



z=0 



0, (36) 



where we have set J^ xt = 0. 

Imposing these boundary conditions we have solved numerically Eqs. (32) by using the COM- 
SOL package [38]. Solutions of this form correspond to straight vortex lines centered at r = 0. 
For n = the only solution we find is the homogeneous solution of Section 4 with B = 0: this is 
the famous Meissner effect of superconductors. In Fig. 5 we give the condensate and the magnetic 
field as a function of r for the n = 1 vortex. We have checked that our solutions respect the 
model-independent properties given in Table 1. 

Then we have studied the prediction of our model for the penetration depth A' and we present 
it in Fig. 6. This constant characterizes, among other things, the behavior of a v far away from 
the center of the vortex (see Table 1) and this allows us to extract it from our solutions. Let us 
take Za{4>) = e 7< ^ for the sake of definiteness. We observe that when 7 is negative enough A' has 
a qualitatively different large \i behavior with respect to the S-BH phase (and to less negative 
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2 4 6 8 10 12 5 10 15 

r r 

L L 

Figure 5: The condensate (left plots) and the magnetic field (right plots) for the n = 1 vortex. Here we 
set v = 1.2, T ~ 0.15/L, Z^cji) = e 7 ^, Z^((j)) = e a< ^ and the kinetic boundary term in (29) equal to zero, 
g/e b -> 0. 

or positive values of 7). Since increasing fiL for a fixed value of TL corresponds to exploring 
the region of small T/T c , we conclude that the superconductors in the dilaton-BH phase have a 
qualitatively different low temperature penetration depth. In particular A' does not increase its 
value more and more as T/T c is decreased, contrary to what happens in the S-BH phase [18]. 
In the low /i region we eventually find that A' diverges at fi = fi c : this is true in any model of 
superconductivity because when fi ~ fi c the condensate becomes small and the Ginzburg-Landau 
theory can be applied to predict that divergence. 

Moreover, we have extracted the prediction for the first critical field H c i, the minimum value 
of the external magnetic field H for which the system is in the vortex phase. This quantity can 
be computed with the text book formula given in Table 1. There en is the U(l) gauge constant of 
the CFT in the normal phase, (O) = 0, normalized in a way that the kinetic term of a stationary 
vector potential in the free energy is J (PxJ 7 ^/ (4el); taking J 7 ^ to be constant we obtain 

1 1 f Zh 

- = - dzZ A (^). (37) 
e o 9 Jo 

We provide H c i as a function of fi in the left plot of Fig. 7. In the high /i//i c (low T/T c ) region 
we find, like for the penetration depth, substantial qualitative differences compared to the S-BH 
case if the exponent 7 is negative enough: while for the S-BH (and for less negative or positive 
values of 7) H cl — > [18], in the other cases H cl remains sizable even at low temperatures. The 
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Figure 6: The penetration depth as a function of fiL. Here we set v = 1.2, T ~ 0.010/L, Za{<P) = e 1 ^ , 
Zip{4>) = 1 and g/eb — > 0. We choose 7 = —4.4 (for which fi c ~ 0.17/L) on the left and 7 = 0.4 (for 
which /i c ~ 0.37 /L) on the right. 




1 2 3 4 5 6 0.2 0.4 0.6 0.8 1.0 

fIL /iL 



Figure 7: Left plot: H c \ as a function of [iL. Right plot: comparison between H C 2 and H c ±. Here we set 
v = 1.2, T ~ 0.010/L, Za(0) = ei^, 7 = -4.4, Z^(0) = 1 and 5 /e 6 -> 0. 

reason why this occurs is because if c i oc eg, whose low temperature behavior depends strongly 
on whether we are in the scale invariant phase, <f> = 0, or in the phase in which scale invariance 
is broken, <fi 7^ 0. Eq. (37) tells us that in the scale invariant phase eo — > and so does H c i, 
but when scale invariance is broken cq remains finite if 7 is negative enough. From Eq. (6) we 
see that the maximum value of 7 for which this is true is —a/2/ [y 2 — 1). Interestingly, these 
negative values of 7 lead to an insulator at low temperatures in the normal phase. It is possible to 
reformulate this discussion in more general terms than exponential (</>): what actually matters 
here is whether the quantity ZjS$)\z=z h is small enough for small temperatures; comparing with 
Eq. (19), we see that when this is the case the normal phase tends to be insulating. In the region 
where \x ~ fi c we have that iJ cl decreases and eventually vanishes at \x = fi c : this is a model- 
independent feature of superconductivity because at the critical point both the superconductive 
and the vortex configurations go to the normal one. 

Finally, we have analyzed the second critical field H c2 , defined as the maximum value of H for 
which we are in the vortex phase. Since \& goes to zero continuously as H — > H c2 , the backreaction 
of ^ on B is so small for H ~ H c2 that one can take H = B and neglect the dynamics of the 
magnetic field when computing H c2 . This means that we can set a Dirichlet condition for A\ at 
the AdS-boundary and that H c2 is the same for superconductors and superfluids. We will describe 
explicitly how we computed this quantity in Section 6, but here we present the final result. In Fig. 
7 we give H c2 and we compare it with H c \. We observe that H c2 is always bigger than H c \ and 
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arg(O) 



vortex dp-behavior 



quantization of B in vortices 



vortex energy 



first critical field 
second critical field 



Superfluid (SF) 



SF current density 

SF velocity potential in the lab frame 



external velocity in the lab frame 



is frozen 



No 



large R 



B cl -"~~2\n(R/U)/R 2R ^°°0 

H C 2 



GL 



Superconductor (SC) 



EM current density 
condensate's phase 



EM vector potential 



large r 



n + a\^fre 



-r/\> 



yes: J drrB 



n 



finite as R — > oo 
H cl = el{F x - F )/2n 

H C 2 = -^FT- 



Table 1: We compare superconductors and superfluids. X' is the penetration depth, a model- dependent 
constant, which is generically different from A (the inverse mass of ai). a\ is another model- dependent 
constant. F n is the free energy of a vortex with winding number n, while £ GL is the coherence length of 
the Ginzburg- Landau theory. We denote with H the external magnetic field and we normalize it such 
that it coincides with the total magnetic field B at = 0. For superfluids we have H = B. See Ref. [18] 
for a derivation of the model-independent properties in this table using effective theory methods. 



therefore there is a finite range of H, that is H c \ < H < for which the superconductor is in 
the vortex phase (the superconductor is of Type II). The holographic superconductor in dilaton- 
gravity shares this feature with the other holographic superconductors studied so far, based on 
pure Einstein gravity, both in ungapped [18] and gapped [22] phases, and therefore it seems, at 
least to some extent, a universal property of holographic superconductivity. 

Since we are dealing with a Type II superconductor we can conclude that for H > H c i, but 
still close enough to H&, the energetically favorable configuration is the single vortex with n = 1, 
while if H increases further, more and more vortices appear until they form a triangular lattice 
at H ~ H c2 



6 Holographic superfluids in dilaton-gravity 

As we have already said, in the limit in which a M is nondynamical a superconductor is equivalent 
to a superfluid. In this section we focus on such limit. In order not to duplicate the symbols, 
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we use the superconductor notation here, but in Table 1 we give the dictionary to interpret the 
various quantities in terms of superfluid physics. In particular we are interested in the case in 
which an external velocity (denoted by aj in the lab frame) is performed on the system. In the 
formulation of superfluidity we are considering all vectors are expressed in the frame where the 
external velocity is zero, unless otherwise stated. 

Since superfluids are characterized by being irrotational we focus on constant external angular 
velocities, where the only nonvanishing component of a« is 

a„ = |r 2 . (38) 

The external angular velocity coincides with B/2. When B is small enough the system is in the 
superfluid state that is irrotational, so the fluid angular velocity is -a v /r 2 and (J v ) ^ 0. At high 
B superfluidity is eventually broken and (J v ) = (normal phase). For intermediate values of B 
the system is in the vortex phase. To holographically describe this situation we therefore consider 
once again the ansatz in (31). 




Figure 8: The condensate (left plot) and the current (right plot) for the n = 1 superfluid vortex. In this 
plot we set v = 1.2, T ~ 0.15/L, B ~ 0, Z A {<j)) = e 7< ^, Z^<f>) = e a(f > . 

Eq. (38) can be viewed as a Dirichlet boundary condition at z — 0. At large 5 values of r, 
r = R, we impose instead 

d r ^P = , cUo = , A v = ^R 2 . (39) 

5 Since the superfluid vortex free energy diverges in the infinite volume limit (see Table 1), we take a large but 
finite-size system: < r < R. 
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The first two conditions correspond to the physical requirement that the field configuration goes 
to the homogeneous superfluid state at large r, while the last one is a simple possibility compatible 
with the form of a v given in Eq. (38). 

We have numerically [38] solved Eqs. (32) with the boundary conditions in Eqs. (33), (34), 
(22), (38) and (39). In Fig. 8 we present the corresponding vortex profiles. 

The minimum value of B at which the vortex starts being favorable, called B c i, is known for 
large R and is given in Table 1. That formula is model independent and so must be true here in 
particular. In the infinite volume limit, R — > oo, we have that B c i — > 0; superfluids are deep Type 
II superconductors and as such they have the properties mentioned at the end of Section 5. 

As we have already stated, the maximum value of B for which the system is in the vortex 
state, B C 2 = H c2 , coincides with that for superconductors. We have computed H c2 as the value of 
B at which ^ — > and provide it as a function of fi on the right plot of Fig. 7. 

7 Outlook 

We have studied superconductors and superfluids through the AdS/CFT correspondence with a 
dilaton in the gravitational theory. Our results are summarized in the introductory section. Let 
us provide some outlook of our work here. 

We have focused on the case in which the charged scalar \1/ is massless, this corresponds in the 
CFT side to an operator of dimension 3. Giving \1/ a mass would change the dimension of O. It 
may be interesting to do so as this would introduce a third function of the dilaton that represents 
a running mass. Along these lines one could then analyze a more general theory with the given 
field content having several running constants, which may be useful to identify model-independent 
features. 

We found that the holographic superconductors studied here are of Type II, like those without 
the dilaton [18, 22] proposed in Refs. [12, 20]. It would be valuable to know if there is a general 
reason for this property or, on the other hand, it is possible to engineer a Type I superconductor. 
A more general theory, similar to that mentioned in the paragraph above, could be a good starting 
point to address this question. 

The (non ideal) insulator/superconductor transition that can be realized between the dilaton- 
BH and the corresponding condensed phase is an alternative to that proposed in Ref. [20], which 
exploits a well known deconfined/confined transition [40, 41]. It is then natural to ask if the 
properties of the latter scenario (see for example Refs. [42, 43, 22, 44]) are shared by the dilaton 
theories discussed here, for the setup corresponding to an insulator /superconductor transition. 

We leave all these questions for future work. 
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